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■ In this paper we revisit the dynamical dark energy model building based on single scalar field 
involving higher derivative terms. By imposing a degenerate condition on the higher derivatives in 

. curved spacetime, one can select the models which are free from the ghost mode and the equation 

■ of state is able to cross the cosmological constant boundary smoothly, dynamically violate the null 
energy condition. Generally the Lagrangian of this type of dark energy models depends on the 

, second derivatives linearly. It behaves like an imperfect fluid, thus its cosmological perturbation 

theory needs to be generalized. We also study such a model with explicit form of degenerate 
Lagrangian and show that its equation of state may cross —1 without any instability. 



PACS nuniber(s): 98.80. Cq. 



I. INTRODUCTION 

Oi 

(— I I The recent data from type la Supernovae (SNIa) and cosmic microwave background (CMB) radiation and so on 
I— I. have provided strong evidences for a spatially flat and accelerated expanding universe. In the context of Friedmann- 
Robertson- Walker (FRW) cosmology with Einstein gravity, this acceleration is attributed to the domination of a 
, component with negative pressure, called dark energy. So far, the nature of dark energy remains a mystery. The- 
' oretically, the simplest candidate for such a component is a small positive cosmological constant, but it suffers the 
[ difficulties associated with the fine tuning and the coincidence problems. Therefore, many physicists are attracted by 
p\| ■ the idea of dynamical dark energy models, such as quintessence [243, phantom k-essence quintom 6] and so 
Tij" I on (see Refs. 0-[Ill for recent reviews). 

. Although the recent fits to the data [l^l , in combination of the 7-year WMAP [ij] , the Sloan Digital Sky Survey 
[3, and the recently released Union2 SNIa data [l^, show remarkably the consistence of the cosmological constant, 
it is worth noting that a class of dynamical models with the equation-of-state (EoS) across —1, dubbed as quintom 
which dynamically violates the null energy condition (NEC), is mildly favored [6|. However, it was noticed that 
consistent single field models realizing the quintom scenario are difficult to be constructed. For example, the EoS of 
■ quintessence is limited to be in the region — 1 < w < 1, while for phantom w is always smaller than —1. It was proved 
/\ ■ that in models described by a single perfect ffuid or a single scalar field with a Lagrangian of k-essence form [1^| , the 
cosmological perturbations encounter a divergence when the background EoS crosses This statement 

was explicitly proven in Ref. [2l| as a "No-Go" theorem for dynamical dark energy models. 

To realize a viable quintom scenario, one usually needs to add more degrees of freedom into the dark energy budget. 
The simplest and also the first quintom model was constructed by a combination of a canonical scalar and a phantom 
scalar field f6j. However theorists are still interested in pursuing single field quintom models. The first single scalar 
field quintom model was realized by introducing higher order derivative terms [22j. also see (23[for generalization. It 
was also obtained in the frame of nonlocal string theory p3 | and decaying tachyonic branes[25| . However, there is 
a quantum instability due to an unbounded vacuum state in such type of models [1^ [13] . A possible approach to 
stable violations of NEC is the ghost condensation of Ref. [l^ , in which the negative kinetic modes are bounded via 
a spontaneous Lorentz symmetry breaking, although it might allow for superluminal propagation of information in 
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some cases |29j . Various theoretical realizations of quintom scenarios and their implications for early universe physics 
were reviewed in Ref. (see also [30|])- 

Recently, a scalar field model which stably violates the NEC has been studied extensively, which is the so-called 
Galileon [31j. It was viewed as a local infrared modification of General Relativity, generalizing an effective field 
description of the DGP model [33] • The key feature of these models is that they contain higher order derivative 
terms in the action while the equation of motion remains second-order in order to avoid the appearance of ghost 
modes, realizing the idea pioneered by Horndeski thirty years ago 33]. Later on, various phenomenological studies 
of this type of models were performed, namely, see Refs. 34-45]. Motivated by the feature of the Galileon model, 
in this paper we revisit quintom dark energy models containing higher derivative terms. We start from a general 
covariant Lagrangian of single scalar field involving higher derivative terms and pursue how to keep the model free 
from extra degree of freedom. We show that it is able to eliminate the ghost mode by imposing a degenerate condition 
that the Lagrangian only depends on the second derivative terms linearly. Based on the degenerate Lagrangian, we 
build an explicit dark energy model and study its dynamics of its homogeneous background and its perturbations. 
Our numerical calculations show that the EoS is able to cross the cosmological constant boundary smoothly and the 
perturbation modes are well controlled when the crossing takes place. We understand the reason of realizing the 
single field quintom without ghosts is that, such a single scalar field model is no longer be able to correspond to a 
perfect fluid. Thus, this model does not conflict with the "No-Go" theorem for quintom dark energy model building 
as proposed in [21|]. 

This paper is organized as follows. In Section II, we simply review the difficulty of constructing a single field dark 
energy model which gives rise to quintom scenario. In Section III, we start with a single field action involving higher 
derivative terms, and present the general analysis. We also discuss under what condition the higher derivative terms 
do not bring a ghost mode to the effective field description in Section IV. Section V is devoted to the study of the 
background and perturbation dynamics of single field dark energy model with degenerate higher derivatives in the 
frame of flat FRW universe. Specifically, we present an explicit model of quintom dark energy with degenerate higher 
derivatives. We perform numerical computation to illustrate such a model can realize the EoS across —1 smoothly. 
Section VI is the summary. 

II. THE DIFFICULTY OF SINGLE FIELD DARK ENERGY MODELS WITH EOS CROSSING 1 

We begin by briefly reviewing the difficulty of constructing dynamical dark energy models with the EoS across the 
cosmological constant boundary. As demonstrated by several groups 0, [T7l - l2l| , the EoS of dark energy based on 
single perfect fluid or single k-essence scalar field (quintessence and phantom are special cases of k-essence) cannot 
cross —I remaining finite perturbations. The key point to the proof of this no-go theorem is that in both scenarios 
the pressure perturbation has a gauge invariant relation with energy and momentum density perturbations. 



where the dot represents the derivative with respect to time, and H = a/a \s the Hubble rate with a the scale factor. 
The sound speed square is defined as the ratio of 5p/5p at the frame comoving with the dark energy, it should be 
positive definite to guarantee the Jeans stability of perturbations at small scales. The momentum density perturbation 
is defined asV = ik''T^ in Fourier space, here is the — i component of the energy momentum tensor. For perfect 
fluid the pressure is a function of the energy density only, p — p{p), and cj. — p/ p, hence Sp ~ {p/p)6p. For k-essence 
field the sound speed square is generally different from p/ p, however, from Eq. ([T]) we can see that in both cases the 
pressure perturbation diverges at the crossing point p = with finite V and density perturbation. As we learn from 
the gravitational field equation, a divergent pressure perturbation will lead to arbitrarily large metric perturbations. 
This instability can also be seen from the equation of motion or the action for the perturbations. A generic form of 
Lagrangian for a k-essence field is only a function of the scalar (p and its first derivatives X — ^(t>V^4), of which 
the action is expressed as. 



where g is the negative determinant of the metric tensor g^^, . Its energy momentum tensor has the same form as that 
of perfect fiuid. 



5p^clSp + 3H{cl-i^)^ , 



P 



(1) 




(2) 



T^-^ = -pg^" + {p + p)ut'u' 



(3) 
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where 

p = C, p = -p + 2Xpx , u^' = . (4) 



We have used the subscript X to indicate the partial derivative with respect to X. The evolution of the k-essence 
perturbation tt — (^(x^ t) — (j)(t) is governed by the equation which may be obtained via variational principle from the 
following action of second order, 



- / (rxdta''[pxTT ^d^TrdiTT +p007r tt J . (5) 

Z J a a"^ 



We have neglected the metric perturbations for simplicity. We can see from the action that the sound speed square 
is = px/px- However, from the relation p + p = 2Xpx we know that px vanishes at the point of crossing and 
changes the sign after the crossing. To guarantee > 0, px should also vanish at the crossing point and change the 
sign afterwards, just as px- The vanishment of px will make the equation of perturbation singular and the amplitude 
of the perturbation arbitrarily large around the crossing point. 

Such kind of instability is classical. Another difficulty emerges when the quantum effects are considered. At the 
phantom phase w < —1, px is negative and > requires px < 0. The action ([5]) showed that the kinetic term of 
TT has a wrong sign. This means tt in this phase is a ghost which brings the problem of vacuum instability due to the 
existence of negative energy states or of violating unitarity by negative norm states. 

III. SINGLE SCALAR DARK ENERGY MODEL WITH HIGHER DERIVATIVES 

To avoid the problem of singular perturbation possessed by single k-essence field, earlier quintom model build- 
ings introduced multi-degree of freedom explicitly. For example in Ref. @ , the quintom model is constructed by a 
quintessence field and a phantom field. Though the total EoS of these two fields crosses — 1 during the evolution, each 
component does not cross this boundary and has regular perturbation. Besides the multi-fluid or multi-field models, 
it is still interesting and important to pursue quintom models with single degree of freedom. To this end, some 
extensions beyond the perfect fluid and k-essence field are proposed in the literature including the higher derivative 
field theory [23 |. the non-minimal coupling to the gravity [46 1, the constrained scalar field which violates Lorentz 
invariance locally [47|, and so on. In this paper we only consider the first extension. 

For a scalar field with higher (but finite) derivatives, its Lagrangian generally has the form, 

£ = £(0,0^1, , (6) 

where = V^^t/), 4>fj.ip,2 = ^p.2^p.i4' ^-nd so on are the covariant derivatives of (f> and N > 2. The equation of 
motion from this Lagrangian is 

f + D-irv„...v.(^) = o. (7, 

Generally this is a 2A'^th order derivative equation, the whole system contains iV degrees of freedom and some of them 
are ghosts. In order to keep the discussions simple and without loss of general properties of higher derivative field 
theories, we only consider the case = 2 in curved spacetime, the Lagrangian is a scalar function of 0, 0^ and (f>p,u- 
The equation of motion is 

dC „ , dC . ^ „ , dC . ^ 
Expanding this equation and considering the symmetry (p^^, = 4>^^^ we have the following equation. 



dC d'^C ^ ^ d'^L d'^C ^ d^L +2—^-^—0 , ^ ^ 

dcj) d(j)d4)f_, ^ d(j)d4)f^u d(j)^d(j)fj^ dcjidcjidcj),^^ ^ " d(j)d(j}pd4>f^^ ^'^ " d4>pd(j)ad4>f^^ 

_d^£_ _ d^C ^ d^C ^ d^C ^ 

777 TTT [vi^PfJ, ~ 'Pi^fJ.p) + 2 „ , „ , ^77 (Ppap.Vi' + ^777 777 777 (Ppap.Vav + 777 777 777 (PajS pV pcru + 
0(t>pd(l)py d(l)d(l)pad4>pu d(j)ad4>pad(l)f,ty d4>apd(l}pad(j)f,i, 

d^C , 



0. (9) 
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If the matrix — is non-degenerate, this is a fourth order differential equation. To solve this equation we need 

to impose the values of 0, 0^, 0^1/ , (j)^ivp at the initial surface. This means such a model essentially possesses two 
dynamical components. 

The first dark energy model of single scalar field with higher derivatives was proposed in (22] where a simple example 
is taken to illustrate the feature of crossing the boundary of —1. The sample model has the Lagrangian of the type, 

L = AX + B{Ud^f -V{^) , (10) 

where □ = V^V, V{(j)) is the potential and A, B are constants. In Ref. [l^l, it was shown explicitly that this model 
is equivalent to the double-field model with lower derivatives, but one field has negative kinetic term. This can be 
shown as follows. Taking an integration by parts the Lagrangian (fTO|) can be rewritten as a function which does not 
depend on X. So it belongs to the more general class in which the Lagrangian is an arbitrary scalar function of (p 
and □(/), 

C = C{<P, Ucj,) . (11) 

Now we use this Lagrangian to illustrate how is the non-degenerate higher derivative field model equivalent to multi- 
field model with lower derivatives. Non-degeneracy means Cbb ^ 0, where the subscript B represents the derivative 
with respect to □</). We first introduce an auxiliary field defined as 

'^ = I^B- (12) 

Due to the non-degeneracy, we may convert the above equation to get Dc^ as a function of </) and Then change the 
independent variables (□(/), ^) to (^i), (p), we accordingly have the Legendre transform of the Lagrangian in Eq. (fTTj) . 

t/(0, = 0n0((/., , (13) 

it may be considered as a potential of (p and (p. After finding the potential [/, the higher derivative term in the 
Lagrangian (fTTj) can be removed with the price of introducing extra field, 

£ = (pn<^ -[/((/,, ^) , (14) 

which is equivalent to the following Lagrangian 

£' = -V^(^V^(^ - C/(0, ^) . (15) 
This equivalent Lagrangian has a more clear form of double fields, 

= iv^0iv^<^i - iv^(/.2V^</)2 - c/(^i, 4>2) , (16) 



through the field redefinitions 



-^(0-(^), 02 = -L(0 + ^). (17) 



The mode 02 is a ghost which violates the null energy condition because its kinetic term has a wrong sign. This 
explains why the model (1111) may cross the cosmological constant boundary without divergent perturbation. However 
at the quantum level, the non-degenerate higher derivative model is plagued by the existence of ghost mode. 



IV. DEGENERATE HIGHER DERIVATIVE MODEL 

In the degenerate higher derivative model, the matrix — in Eq. Q is identically zero. Correspondingly, all 

the third and fourth order derivative terms in the equation of motion (|9]) disappear, i.e.. 
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where we have used the commutation of the covariant derivatives 4>upp, — 4>v^i,p = 4'(tR'^ i^^p with R" the Riemann 
tensor. The equation of motion remains to be a second order differential equation, but there is a curvature-field 
coupling term appeared in it even though we only consider the minimal coupling to the gravity in the Lagrangian. 
This means the degenerate model has no extra degree of freedom. ^ 

Zero matrix „ , ^ — implies the Lagrangian only depends on the second derivative terms linearly. With Lorcntz 
invariance, £ can only be 

C = K{(l), X) + G{(t), X)U4) + X)Vt'(tN''(tN^,V,,(t> . (19) 

The box term at the right hand side with G = X was considered in the context of Galileon theory 'si'] and its 
generalization was studied in [ssj , named as KGB model and in inflation model building 43] , named as G- inflation. 
The box term is equivalent to --2XG^ — G x'^ ^ 4>^'^ 4>^ p^ iy<t) after integration by parts and dropping a surface term. 
By redefinitions of K{4), X) and F{(t) ,X), the degenerate Lagrangian may be generally written as 

C = K{(t), X) + F{(t), X)V'"0V''0V^V^0 + i^V^XV^0 . (20) 

Now we will investigate whether the dark energy model from this Lagrangian can stably cross the boundary of 
cosmological constant by studying its background evolution and properties of perturbations. 
With the notations (|20p the equation of motion becomes 

- 2XKx^ - KxO^ - XxxV^XV^.^ + F[{a(j)f ~ V^V,0V''V"<^ + i?^,V'>V>] 
+ 2F^(V^XV^</. + 2XD<p) + 2Xi^^xV^XV^0 + 4X^F^^ + i^x(D0V^XV^</. - W^.XWX) = . (21) 

The energy momentum tensor which sources the gravitational field is obtained through the variation of the action 
with respect to the metric tensor, 

^,.^_2S^^ (22) 

for the degenerate model it is 

T^'^ = -{K + FVpXVP<j))g'"' + {Kx - 2XF^ - FU(t))yf'(tN''(t) + F(V^XV''0 + VXV^.?!)) . (23) 

We can read off the pressure and energy density from this energy momentum tensor up to linear order of perturbations 
around the homogeneous background, 

p = K + FV^XV>'(j) , (24) 
p = -K + FVpXV>'(l) + 2X{Kx-2XF^-FU(j)) . (25) 

In terms of the fluid variables p and p, the energy momentum tensor may be rewritten as 

T^'' = -pg^'^ + {p+ p)ut'u'' + (2X)3/2F(a^u'^ + a^u^") , (26) 

where by analogy with k-inflation [4^ or k-essence Q we have defined the four velocity = V^4>l\j2X which is 
normalized as = \,a^ = pu'^ is the four acceleration which is orthogonal to the velocity, i.e., a^u^^ = 0. We 

have also used the relation 1361: 



V^X = 2Xa'^ + uPVpXu'' . (27) 

This energy momentum tensor (j26p does not have the form of perfect fluid due to the last two terms depending on the 
four acceleration. In the language of relativistic imperfect fluid, can be identified as heat flow. Further, one can 



^ In the flat spacetime it is not necessary to require — = to keep the equation of motion at the second order, for example if 

dtb^ 8d> — ~ C {rj'^'^ ri'"' — rj^Prj"") with constant C in the Minkowski space, all the third and fourth order derivative terms vanish in the 



equation of motion because (pf^^pa- is totally symmetric under the interchanges of the indices. See Ref. 38] for some more discussions 
about the case in flat spacetin 
This has be shown in Ref. [33 



about the case in flat spacetime. But in the curved spacetime, „ ^ ^ — = is the unique way to discriminate higher derivative terms. 
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prove straitforwardly that a'' is space like and its zero-th component vanishes at both background and Unear levels, 
its spatial components should be first order variables. So if we redefine the four velocities as 

= u'' + ^ '- , (28) 

P + P 

the energy-momentum tensor should be 

T^"" = -pg^""" + {p + p)u^'u'' - I^^l-^a^'a'^ ~ -pg^''' + {p + p)u^'u'' . (29) 

P + P 

In the last step, we have considered the fact that the product a^a'^ is a second order variable and can be neglected if we 
restrict our studies on the background evolution and the linear perturbation theory. This equation means the energy 
momentum tensor of the degenerate model has apparently the same form of perfect fluid if we neglect perturbations 
of higher order. But it is essentially different from perfect fluid, especially the relation ([Ij between the pressure and 
density perturbations is lost in this model. This is the very reason why the degenerate higher derivative model is 
possible to realize the quintom idea and avoid the problems possessed by single k-essence field. 



V. COSMOLOGY WITH DEGENERATE HIGHER DERIVATIVE DARK ENERGY MODEL 

In this section, we consider in more detail the dynamics of the dark energy model (j20p . The background universe 
is a spatially flat FRW spacetime, in which the metric is given by 

ds^ = dt^ - a^{t)S,jdx'dx^ . (30) 

At the background level the field (j) is homogeneous, so we soon have 

p= -K + 2X{Kx - 2XF^ - 3iJ0F) 
p = K + 2X4)F 

p+p = 2X{Kx - 2XF^ ~ m^F + '<i,F) , (31) 

where we have considered X = 0^/2 , V^XV^0 = 2Xip and □</) — (f) + 3H<j). The energy conservation law in the 
expanding universe p + SH{p + p) = is identical to the equation of motion (|2ip . For the velocity u'^ only the time 
component is non-zero, = 1, — 0. The key point for this model to cross the cosmological constant boundary is 
that Kx — 2XF^ — 3H(/)F + (pF could evolve from the positive region to the negative region or vise versa providing 
the energy density always positive. In addition, we have to check whether the perturbations are stable. 

Similar to the analysis in single k-essence model, the perturbation of the scalar field is a small deviation from the 
homogenous background 

(j){x, t) = 0(i) + 7r(f, t) . (32) 

For complete consideration wc should also include the perturbations of spacetime. However, for dark energy, it is 
believed to be subdominant in the universe for most time and had tiny contribution to the curvature. So for studying 
the dark energy perturbations, it is safely to neglect the metric perturbations which are sourced mainly by other 
matter. This approximation will greatly simplify the analysis. There are two ways to get the linear perturbation 
equations. One is to directly perturb the equation of motion (PT|) around the background evolution. Another one is 
to adapt the variational principle from the action which is second order of perturbations as we have shown in Sec. II 
for k-essence field. Here we will use the second way to discuss the perturbation of the scalar field. For this purpose 
we firstly expand the action 

S = j d'^Xy/g{K + FV^Xy^'cj)) (33) 

to the second order of the perturbations. After straightforward but tedious calculations, we obtain the desired second 
order action for tt, 

S'(2)(7r) ^ I [ d^xdta^iAn^ - ^d.nd.TT + Ctt^) , (34) 
2 J 
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with 

A ^ Kx + 2XKxx - GH^iF + XFx) ~ &XF^ - 'iX^Fx^ , 
B = Kx-2{^ + 2H4>)F - AXF^ - 2X'<}>Fx , 

C = j^{6HXF^ + 2X^F^^ ~ ^Kx4) + iH{&HXF^ + 2X^F^^ - 4>Kx^) + 2X'^F^^ + K^^ . (35) 
We can see from the action that the sound speed square is defined as 

^2^B^ Kx - 2i4> + 2H^)F - AXF^ ~ 2X'^Fx 

A Kx + 2XKxx-GH^{F + XFx)-8XF4,-4X^Fx4>' 

If the universe is dominated by the scalar field, as discussed in the KGB model [35] or the G- inflation model [4^], a 
full treatment of the gravity-(/) coupled system based on the (Arnowitt-Deser-Misner) ADM method will give a slightly 
different sound speed squared, 

^,2 _ Kx - 2{4> + 2H^)F - AXF^ ~ 2X'4>Fx ~ 2X^FyM^ 

^ Kx+ 2XKxx - 6H4>{F + XFx) - SXF^ - AX^Fx^, + GX^F^M^ ' 

where both numerator and denominator are modified by terms suppressed by the Planck mass — I/SttG. This 
sound speed depends on the gravity theory, here the gravity theory is Einstein's general relativity. For dark en- 
ergy studied in this paper we will consider the sound speed in Eq. (1361) to express the propagating velocity of the 
perturbations. 

The classical stability requires > 0. Furthermore, the absence of ghost mode corresponds to 

A = Kx + 2XKxx - 6H^{F + XFx) - 8XF^ - 4X^Fx0 > . (38) 

Compared to Eq. pil) one can find that neither A nor B is proportional to p + p, hence when the dark energy crosses 
the cosmological constant boundary, p + p = 0, both coefficients A and B are not vanished in general. The equation 
of motion would be regular at the crossing point. This is different from the case of k-essence model discussed in 
Section II. So in this single field model, for particular choices of the functions K and F and corresponding model 
parameters, it is possible to find solutions in which the equation of state of dark energy evolves across —1 but both 
and A remain finite and positive as we will show explicitly below. Within these solutions the fluctuation of the scalar 
field has the right kinetic term to circumvent the pathology of ghost, even though its background part has w < —1 
violating the null energy condition. 

In order to illustrate the realization of quintom scenario explicitly, we study an explicit form of the degenerate dark 
energy model. The Lagrangian we are considering is simple: 

C = -X- ciV^XV^'cj) + C2X(j)'^ , (39) 

where ci^2 are constants. Compared with the notations in Eq. (j20l) . one may find that K = —X + C2X(j)'^ and 
F = — ci. Note that the third term can also be viewed as an "effective mass term". From Eqs. ([M]) and we get 
the pressure and energy density of this model respectively as: 

p = {c2(f>^ -l + 6ciH4>)X , (40) 
p - {C24>^ - 1 - 2ci'4>)X . (41) 

At the point where the equation of state crossing —1, we have p + p = Q, i.e., 

ci(<ii-3i?(/.)=C2(/'2~l . (42) 

In addition, we have the equation of motion for the scalar field, which comes from Eq. (1211) . 

(1 - QciH^ - C2(t)^)4> + 3H{1 - C2(t)^)^ - (c20 + 9ciij2 + 3ciH)^^ = . (43) 

We performed the numerical calculations of the background evolution and the coefficients for perturbations for 
different parameter choices in Figs. [T]and[21 In both figures, we have used the unit SttG = 1 and chose the parameter 
ci = 3.0 X 10^^^. The choices of parameter C2 are different for these two cases. In the first figure C2 = and in the 
second one C2 = 0.05. In both cases we set the initial conditions as = —1.2 and (pi = 10~^^ well within the matter 
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dominated era in order to make the evolutions of dark energy consistent with the observations. We can see that in 
both cases the equation of state of dark energy crosses the cosmological constant boundary, and the present values 
of wde and ^Ide fits well with the observational data. At present time Ina = 0, w^e — —1.158 and —1.167 for the 
first and second cases respectively. These are consistent with the current constraint by the observational data from 
CMB, LSS and Supernovae. For example the analysis in Ref. ^49^] gives wde = —1.143 ± 0.160. The values ofil.DE m 
both cases are around 0.725, also compatible with the observations [31 . Moreover, the parameter A defined in ([35]) 
and the sound speed squared defined in (|36p are both positive during evolutions in these two cases. This means that 
at the crossing there is no ghost instability and the perturbation is regular and classically stable. 

Finally we should emphasize that with different parameters, these two cases give rise to predictions of different fate 
of dark energy in the future. For the first case where C2 = which reduces to the case discussed in [ll], wde crosses 
— 1 and approaches —2 at early time, then when the dark energy dominates the universe wue will approach to — 1 
forever and the universe enters into the De Sitter phase. This can be understood as follows. With vanished C2, the 
model is symmetric under the field shift —?>(/) + C with constant C. The corresponding Noether current is 

= ciV^X + V''(/)-cin0V'^(^ , (44) 

and the conservation law = is nothing but the equation of motion. For the background, cf) is homogeneous 

the conservation law simplifies as 

j° + 3i? J° , (45) 

and the "charge density" = 4> — SciHip^, which scales as a^'^. From it we can solve for (j> that 

. 1 ± VI - 12ciH J" 
6^ • 

On the other hand the energy density and pressure are respectively 

p = (_i + 6cii/</>)^ 

p = -(l + 2ci^)^ . (47) 
To guarantee the positivity of the energy density, GciHif) > 1, in Eq. pS)) . we should choose the plus sign, i.e.. 



bciH 

^ soon after the beginning J*^ approaches to zero so that (j) ^ 1^777 't' ~^ ~3cTh^- This 
time the universe is still dominated by matter and the Hubble rate H scales as a~^^^, it is easy to obtain that 
approaches to the constant ^ . So the equation of state 

WDE = ^ -(1 + 2C10) = -2 . (49) 

— 1 + 6ciH(p 

This is a phantom attractor solution to the model. At this period, the energy density of the dark energy increases as 
fast as and the matter energy density decreases as a~'^. So with the initial conditions we have chosen the universe 
transits to the phase of dark energy domination at low redshift. In this phase is more close to than in the 

matter dominated era. However, the Hubble rate is 

H^^li-,^,c,Hi,)^^^^^. (50) 

SO H approaches to the constant (54c^)^^/^ and ^ — ?> 0. The equation of state wde ^ ^(1 + 2ci0) — > —1. The 
universe enters into the de Sitter phase in the future. 

For the second case where C2 — 0.05, the shift symmetry is violated by the "effective" mass term. The "charge 
density" satisfies the equation 

i° + 3HJ° = -C20^<^ , (51) 




FIG. 1: The equation of state Wde, the energy density ratio ^Ide, parameter A and sound speed squared of dark energy of 
the model (|39() . The parameters are chosen to be ci = 3.0 x lO^'^^, C2 = and the initial conditions are chosen in the far past 
when the universe was matter dominated. The equation of sate crosses —1 recently, and approach to —1 in the future. Both 
A and are positive during the evolution. The universe will be dark energy dominant in the future. 



where the density J" gets a minor modification compared with the first case, 

j" = (l-C20^)(/)-3cii?(/.2 . (52) 

With the parameters and the initial conditions we have chosen, the effective mass term is very small at the time of 
matter domination. The early evolution of dark energy in this case is similar to the first case, its equation of state 
crosses —1 and approaches to the phantom attractor solution wde = —2. Because the initial value of (p we have 
chosen is negative and (f) ~ ^ciH 't' increases from the negative region to the positive one, so the effective mass 
term remains small during a fairly long period. Recently when the universe shifts to the dark energy dominated 
phase, Wde increases from —2 and approaches —1. However, the value of the field 4> increases continually and the 
effective mass term becomes more and more important. So the universe will not enter into the de Sitter phase and 
Wde increases continually and gets back to larger than zero again. When C2(f>^ increases to the value much larger than 
unity, it will dominate the energy density and the pressure of the dark energy as seen from Eq. (|40p . The equation 
of state Wde — > 1, but this phase is not stable because the energy density of the dark energy will decrease as fast as 
a~^. Hence in the future dark energy might exit domination and the universe will return to matter-dominant stage. 
This phenomenon of the model has not been discussed before. 



VI. SUMMARY 



The model building of dark energy which is possible to cross the cosmological constant boundary, i.e., the quintom 
dark energy has attracted a lot of attention in the literature (for example see Refs. [13] )• A very interesting question 
is how to build a single scalar model which gives rise to the crossing without any instabilities. The past studies have 
shown that a single scalar field satisfying a generic k-essence Lagrangian cannot give rise to the stable violation of NEC. 
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FIG. 2: The equation of state wns, the energy density ratio Qde, parameter A and sound speed squared of dark energy 
of the model (|39p . The parameters are chosen to be ci = 3.0 x lO'^'^'^, C2 = 0.05 and the initial conditions are chosen in the 
far past when the universe was matter dominated. The equation of sate crosses —1 recently, and becomes larger than in the 
future. Both A and ci are positive during the evolution. The universe will be matter-dominant in the future. 



However, it becomes possible when higher derivative operators are introduced. The non-degenerate higher derivative 
model can eliminate the classical instabilities but is still plagued by the ghost mode. In this paper we considered the 
degenerate higher derivative model inspired by the Galileon theory and its generalizations. We started with a single 
scalar of which the Lagrangian contains generic higher derivative operator, and then suggested a degenerate condition 
to eliminate the extra degrees of freedom brought by these higher derivative operators. Using the degenerate condition 
in the curved spacetime, the Lagrangian can be reduced to the form ((20)) . We studied the cosmological application 
of this model, including the background and perturbation analysis. Particularly, we chose an explicit example and 
performed a detailed numerical computation. Our numerical results verified that this model is able to realize the EoS 
to cross the cosmological constant boundary and the perturbation evolves smoothly without any pathologies. 

A quintom model which successfully violates the NEC without leading to quantum instability also has important 
implications to the early universe physics. In Refs. [5l|, the authors found that a quintom model can avoid the big bang 
singularity widely existing in standard and inflationary cosmologies. In this picture, the moment of initial singularity 
can be replaced by a big bounce. Based on this scenario, the correspondin g pe rturbation theory has been extensively 
developed, namely, on adiabatic perturbations [s^ - ls^ ]. non-Gaussianities |55l| . entropy fiuctuationsfssl. |56|. and the 
related preheating phase [57'|. Recently, it was observed that the Galileon model with the EoS across —1 exactly leads 
to a bouncing solution in the frame of the fiat FRW universe[5l|. Consequently, we expect that bouncing cosmologies 
can be realized in a generic quintom model with degenerate higher derivative operators. 
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